Gain narrowing in few-atom systems 
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Using a density matrix approach, we study the simplest systems that display both gain and 
feedback: clusters of 2 to 5 atoms, one of which is pumped. The other atoms supply feedback 
through multiple scattering of light. We show that, if the atoms are in each other's near-field, the 
system exhibits large gain narrowing and spectral mode redistribution. The observed phenomena 
are more pronounced if the feedback is enhanced. Our system is to our knowledge the simplest 
exactly solvable microscopic system which shows the approach to laser oscillation. 
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In a laser, light is generated by a combination of light 
amplification by stimulated emission and optical feed- 
back IjJ. In order to study the basic physics of these 
processes, there has been an intensive search for laser op- 
eration in fundamental systems |2|. The resulting drive 
toward miniaturization has led to, among others, the re- 
alization of vertical-cavity semiconductor lasers j^] , dye 
microsphere lasers j^ , microring and microdisk semicon- 
ductor lasers [5| and photonic bandgap lasers e| . As laser 
systems are made smaller, a purely macroscopic descrip- 
tion becomes inadequate and microscopic considerations 
should be taken into account. An interesting example of 
lasers which require a (partially) microscopic treatment is 
the class of one-atom lasers 7], in which the gain medium 
is reduced to a fundamental level, while macroscopic mir- 
rors provide feedback. Another, contrasting example is 
the class of random lasers 8] , in which optical feedback is 
provided by scattering from microscopic particles, while 
the gain medium remains macroscopic. Obviously, nei- 
ther the feedback mechanism nor the gain medium can 
be reduced to less than one atom. 

In this Letter, we explore the most fundamental sys- 
tem displaying both gain and feedback: a single pumped 
atom, surrounded by one or more passive atoms provid- 
ing optical feedback by scattering. The atoms are posi- 
tioned in free space in the absence of a cavity. The ab- 
sence of a cavity and its modes differentiates our model 
from models in which atoms interact via a single field 
mode such as, e.g., atoms in a single-atom mascr 9']. Our 
few-atom system can be described fully microscopically, 
without any quasi-classical or paraxial approximations. 
We show that this system, though very simple, shows 
surprisingly strong spectral gain narrowing and mode re- 
distribution, indicating an approach to laser oscillation. 
We further demonstrate that the observed phenomena 
are more pronounced as the number of atoms increases, 
in correspondence with the intuitive TV — > cx) limit. 

The "atoms" could be implemented as any type of 
sub-wavelength quantum objects, for example: trapped 



atoms, quantum dots ,10|, trapped ions |ll| or dye 
molecules. Each atom interacts with the electromag- 
netic field by its transition dipole moment, which results 
in scattering of light. One of the atoms is continuously 
pumped, causing it to not only scatter, but also amplify 
the light. This iV-atom system has optical feedback due 
to the fact that in the process of stimulated emission, a 
single atom scatters the stimulated photon isotropically, 
in contrast to the general notion that stimulated emis- 
sion preserves the "direction" of the photon, which is 
only true in a macroscopic gain medium |l2|. 

We assume the atoms to be fixed in space, e.g., by a 
tight trap or a solid matrix. Each atom has three relevant 
energy levels: the ground state a, a highly excited state 
b and the upper state of the relevant c ^> a transition 
c, as depicted in Figure ^ (generalization to a four- level 
system is straightforward) . One of the atoms is pumped 
with light which is resonant with the & — > a transition. 
The pump intensity is expressed by the dimensionless 
parameter W — ri^/FbcTca , where Q, = jdba • Epump\/^ 
is the Rabi frequency of the pump field, with d^a the 
dipole moment of the b ^ a transition |l3| . The decay 
rates from 6 to c and from c to a are given by The and 
Tea respectively. We can consider the atoms to be effec- 
tive two-level (a-c) systems if decay from c to a is much 
slower than from 6 to c and decay from 6 to a is negligible 
compared to other decay processes, as is usually the case 
in a three-level laser. 

The ensemble-averaged populations of the atomic lev- 
els and the coherences between them are described by a 
density matrix 14] (7(i), whose evolution is given by the 
master equation 



d 
di 



a = Cnd^ + -CdCT, 



with as non-dissipative operator 
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FIG. 1: (color online) The three-level system a h c. Fjc and 
Tea are the decay rates from & to c and c to a respectively. 
Decay from 6 to a is negligible compared to other decay pro- 
cesses. The red dashed arrow expresses the interaction with 
the pump field. 



•-<—►• 




and as dissipative operator 



1 ^ 



w 
T 



Tea (crS^-S^+ - S^+aS^-^ + H.c. (3) 



The pumped atom is labeled "1" in expression ©. The 
operators 5"*+ and S"^^ respectively raise and lower the 
state of atom i. The coupling between different atoms is 
quantified by 



1 T c 

2 i^r.a 



Goi^ca,rrn - r„) • /X„, (4) 



for the off-diagonal elements m ^ n. The diagonal el- 
ements of the coupling are given by F"" = Fca,Vn and 
(5"" = LUca = i^c ^ i^aiVn. The tensor Qq represents the 
free-space dyadic Green function [13 , c is the free-space 
speed of light, r^ the position vector of atom i and fii is 
the normalized transition dipole moment of atom i. The 
master equation JQ) is derived by an integration over the 
multimode electromagnetic field to which the atoms cou- 
ple. This integration results in the manifestation of the 
effective coupling term (@J. 

Solving the master equation (^ requires the inversion 
and diagonalization of the associated 2^^ x 2^^ matrix. 
We solve the computer-generated symbolic master equa- 
tion of the total atomic system to calculate the spectral 
distribution of the emitted light [13 . By applying a mas- 
ter equation, we ensure that both elastic and inelastic 
scattering of photons is taken into account in all scatter- 
ing orders. Herein lies one of the major benefits of the 
method we use: the atomic saturation due to inelastically 
scattered photons is significant [12|] and very difficult to 
incorporate in a classical scattering formalism jl^. 

The spectral distribution of the emitted light depends 
on the atoms' spatial configuration and the orientation 
of the transition dipole moments. We focus on configura- 
tions of atoms with interatomic distances L of the order 
c/u>ca- For much larger distances, the feedback provided 
by the passive atoms is limited and only a very small frac- 
tion of the photons emitted by the pumped atom will be 
scattered. On the other hand, for distances much smaller 



FIG. 2: (color) The normalized far-field angle- averaged spec- 
trum versus frequency (in units of Tea)- The inset shows the 
configuration for which the spectrum is evaluated. The pas- 
sive atoms (blue) are positioned in an equilateral triangle, 
with the pumped atom (red) in the center. All transition 
dipole moments are perpendicular to the plane of the atoms. 
The distance L was chosen Q.lc/uica- The spectrum is shown 
for a low pump intensity W — 1.77 (blue) and a high pump 
intensity W = 10.10 (red). 



than the resonance wavelength, the photons emitted by 
the pumped atom will significantly saturate the passive 
atoms and effectively reduce the system's feedback. Con- 
sequently, there is a distance range of the order of the 
resonance wavelength for which the feedback provided 
by the passive atoms is optimal. 

The average atom-photon interaction time is of the or- 
der T~^ , while the time it takes for photons to propagate 
from one atom to another is of the order uj~a ^ ^7a- 
Hence, the information and energy in the system will be 
stored as atomic excitations rather than electromagnetic 
excitations. The system's storage capacity is thus deter- 
mined by the number of atoms. 

The spectral information of the emitted light can be de- 
duced from the Fourier transform of the field-correlation 
function 



gW(T)^(:£;-(t + r).^+(t) 



(2Af)! 



(5) 



in steady-state, where the colons denote normal and time 
ordering for the field operators. As a general result of the 
master equation, the spectrum of the emitted light can 
be expressed as a sum of ■ 



(jv+i)!(jv-i'i! Lorcntzian contri- 
butions [18| . Each of these is characterized by a central 
frequency, a spectral width and a spectral weight. The 
latter expresses to what extent each contribution dom- 
inates the total spectrum. If we increase the pump in- 
tensity, different modes will be subject to different gain 
and, consequently, modes will compete for the available 
population inversion in the system. 

To indicate the effect of an increasing pump inten- 
sity on the spectrum of the emitted light. Figure El shows 
the far-field spectrum for a four-atom configuration. The 
spectrum is averaged over a Att solid angle for a typical 



2ofEl 



Gain narrowing in a few-atom system 



Tom Savels, Allard P. Mosk and Ad Lagendijk 



low {W = 1.77) and a high {W = 10.10) pump intensity. 
At low pump intensity, the emission spectrum is broad, 
while at higher pump intensities we observe a significant 
spectral narrowing. 

If we wish to compare different configurations of atoms, 
we need to visualize the degree of observed gain narrow- 
ing. We proceed along the path originally considered by 
Shawlow and Townes (191] and compare the line width 
of the emitted light to the photon emission rate. We 
therefore determine the full width at half maximum Auj 
of the far-field spectrum I{oj) averaged over a An solid 
angle. We then evaluate the spectral weight /^^ I{ijj)duj 
within the range Aw, which yields the total photon emis- 
sion rate nTca emitted in the range Aoj. The number n 
can then be interpreted as the number of excitations in 
a cavity with decay rate Tea- Since the passive atoms 
can store one excitation each, the average rate nTca de- 
tected in the far-field cannot exceed NTca- Additionally, 
the average photon emission rate is limited by the rate 

WTca\S^~S^^j at which pump photons are absorbed 
by the system. Figure is a parametric plot showing 
the resulting spectral width Acj versus 1/n for the same 
parameters as in Figure [21 for increasing pump intensity. 
We observe a large decrease of Aw, accompanied by an 
increase in n. Since n cannot increase indefinitely, there 
is a critical pump intensity at which saturation of the 
passive atoms sets in, as shown in the inset of Figure |3 
The maximum value of n is relatively low compared to N 
due to the weak coupling between the atoms. Around the 
saturation point, the spectrum broadens while n remains 
locally constant. If the pump intensity increases beyond 
the saturation point of the passive atoms, the emission 
rate decreases. This effect is due to power broadening, 
inherent to the three-level pumping scheme: at the sat- 
uration point of the passive atoms, an increase of the 
pump leads to a decoupling of the pumped atom and the 
passive ones, resulting in the observed decrease of n. 

The observed dependence of 1/n on Aw below satura- 
tion is similar to the behavior found in many macroscopic 
lasers. As was shown by Shawlow and Townes and gen- 
eralized by many others 20], the quantum- limited laser 
line width due to diffusion is inversely proportional to the 
number of photons in the laser mode. It is striking that 
our simple microscopic system exhibits a similar behav- 
ior while we are not in the regime typically considered in 
the Schawlow- Townes relation. 

If the number of atoms increases, the system's storage 
capacity grows. Hence, the maximum value Umax which 
n can attain increases with N . For each given number 
of atoms, there is an infinite number of possible config- 
urations in which the atoms can be positioned. Since 
the dipole-dipole coupling depends on the configuration, 
both rimax and the corresponding width Awmi„ will, for 
a given N , vary with the geometry. 

In order to compare different configurations of atoms. 
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FIG. 3: (color online) The photon emission rate (relative to 
Fca) in the range Alj versus Aoj (in units of Fca). The pump 
intensity W ranges from 1.76 to 13.43, while the inset focusses 
on the behavior around the saturation point. The arrow de- 
notes an increase of the pump. The same configuration as in 
Figure |5| is used, with L = 0.7c/LJca. The red dashed lines in 
the inset denote the position of Aci^mm and n~ 
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we determine how many excitations can be stored with a 
given coherence time Aw^^j^ . For each number of atoms 
A^, we consider those configurations which attain their 
saturation point at a given value of Awmi„ . We then de- 
termine the corresponding number of excitations Umax- 
Figure0]shows the calculated Umax for three different val- 
ues of the coherence time. For every number of atoms, 
different configurations exist which yield the same satu- 
ration value Aujmin- In general, such configurations each 
have a different rimax associated with them, as repre- 
sented by the identically colored symbols. We see that 
the effect of an increase in N is twofold. First, we observe 
that, for a fixed coherence time, the maximum number of 
excitations increases with A''. This trend indicates that, 
as the storage capacity of the system grows, more pho- 
tons with a given coherence time Aw~J„ can be emitted 
by the system. Second, when comparing different val- 
ues of Aujmin in Figure 01 we see that, if the required 
coherence time increases, a larger capacity is needed to 
attain a given number of excitations Umax ■ This relation 
between the number of excitations and the storage ca- 
pacity is in accordance with the intuitive limiting case 

A^ -* OO, l/Umax ^ and AuJrmn ^ 0. 

The efficiency with which incident pump photons are 
converted into photons in the range AtOmin is given by 
the ratio of the output rate Umax^ca and the input rate 
lyFca (S^~ S^^j. The numerical value of the efficiency 
depends on the number of atoms and the configuration, 
but we find as a general trend that the efficiency increases 
with N. For Aw;;^^^ = OASTja considered in Figure JH 
for example, the efficiency increases from 20% for A = 2 
to 24% for A = 5. This indicates that adding more atoms 
leads to a better photon confinement, as we expect. 

The mechanism governing light amplification in the 
few-atom systems presented here is stimulated emission. 
The spontaneous radiation emitted by the gain atom is 
scattered by the passive atoms; subsequent interaction 
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FIG. 4: (color) The relation between the maximum excita- 
tions number Umax and the number of atoms A'". Multiple 
symbols for a fixed A'' represent different configurations. Fixed 
colors represent fixed coherence times. The green triangles are 
for a coherence time 0.47r~;;^. The red diamonds are for a co- 
herence time 0.43r^t,^. The blue circles are for a coherence 
time 0.29r"^ 



with the gain atom generates stimulated emission. The 
presented physical processes correspond to the behavior 
expected of a sub-threshold bad-cavity laser [2^ |23 ■ 

From an experimental point of view, we are convinced 
that our model is relevant for a wide range of realizations, 
such as laser cooled trapped ions |2j|. Selective excita- 
tion of the ions can be achieved by tuning the polarization 
of the pump field and the positioning of the ions. An- 
other possible experimental path would be to implement 
the atoms as quantum dots |22| which form, if bound 
to DNA, bioconjugated superstructures [23 ■ While the 
energy transfer in those structures is somewhat different 
from the one presented here, we anticipate the gain nar- 
rowing phenomenon as presented in this Letter to remain 
conceptually valid. As a third experimental realization, 
we expect our results to stimulate experimental work on 
cold atoms |2J| interacting with pumping fields. Com- 
pared to our current system, the number of atoms in a 
cold-atom cloud is very high; we therefore trust the phe- 
nomena described above to be much more pronounced. 

In conclusion, we described the simplest microscopic 
system which shows both gain and feedback. A sys- 
tem of only a few atoms in each other's near-field shows 
large gain narrowing and mode redistribution. Surpris- 
ingly, our system qualitatively exhibits behavior similar 
to macroscopic lasers. Adding more atoms to the system 
enhances the observed phenomena and allows more pho- 
tons to propagate with a given coherence time. 

Various extensions to our model are possible, among 
which schemes where more than one atom is pumped. 
The challenge of our method lies within the exponential 
scaling of the matrices involved; other approaches such 



as stochastical wave function calculations |2g might help 
overcome the current limits. Conversely, our current sys- 
tem may serve as a building-block for few-body contribu- 
tions in an effective medium approach to, e.g., coherent 
backscattering [l5j. Importantly, we expect our results 
can be tested experimentally, and we are working towards 
a physical realization of our model. 
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